February 2, 2008 17:48 



WSPC - Proceedings Trim Size: 9in x 6in 2007'AMNLWP 



1 



THE D'ALEMBERT-LAGRANGE PRINCIPLE FOR 
GRADIENT THEORIES AND BOUNDARY CONDITIONS 

H. GOUIN 

Universite d'Aix-Marseille, 13397 Marseille Cedex 20, France 
E-mail: henri. gouin@univ-cezanne.fr 

Dedicated to Prof. Antonio M. Greco 



Motions of continuous media presenting singularities are associated with phe- 
nomena involving shocks, interfaces or material surfaces. The equations repre- 
senting evolutions of these media are irregular through geometrical manifolds. 
A unique continuous medium is conceptually simpler than several media with 
surfaces of singularity. To avoid the surfaces of discontinuity in the theory, we 
transform the model by considering a continuous medium taking into account 
more complete internal energies expressed in gradient developments associated 
with the variables of state. Nevertheless, resulting equations of motion are of an 
higher order than those of the classical models: they lead to non-linear models 
associated with more complex integration processes on the mathematical level 
as well as on the numerical point of view. In fact, such models allow a precise 
study of singular zones when they have a non negligible physical thickness. This 
is typically the case for capillarity phenomena in fluids or mixtures of fluids in 
which interfacial zones are transition layers between phases or layers between 
fluids and solid walls. Within the framework of mechanics for continuous me- 
dia, we propose to deal with the functional point of view considering globally 
the equations of the media as well as the boundary conditions associated with 
these equations. For this aim, we revisit the d'Alembert- Lagrange principle of 
virtual works which is able to consider the expressions of the works of forces 
applied to a continuous medium as a linear functional value on a space of test 
functions in the form of virtual displacements. At the end, we analyze exam- 
ples corresponding to capillary fluids. This analysis brings us to numerical or 
asymptotic methods avoiding the difficulties due to singularities in simpler -but 
with singularities- models. 



1. Introduction 

A mechanical problem is generally studied through force interactions be- 
tween masses located in material points: this Newton point of view leads 
together to the statistical mechanics but also to the continuum mechanics. 
The statistical mechanics is mostly precise but is in fact too detailed and 
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in many cases huge calculations crop up. The continuum mechanics is an 
asymptotic notion coming from short range interactions between molecules. 
It follows a loose of information but a more efficient and directly computable 
theory. In the simplest case of continuum mechanics, residual information 
comes through stress tensor like Cauchy tensor 1,2 . The concept of stress 
tensor is so frequently used that it has become as natural as the notion 
of force. Nevertheless, tensor of contact couples can be investigated as in 
Cosserat medium 3 or configuration forces like in Gurtin approach 4 with 
edge interactions of Noll and Virga 5 . Stress tensors and contact forces are 
interrelated notions 6 . 

A fundamental point of view in continuum mechanics is: the Newton sys- 
tem for forces is equivalent to the work of forces is the value of a linear 
functional of displacements. Such a method due to Lagrange is dual of 
the system of forces due to Newton 7 ' 8 and is not issued from a variational 
approach; the minimization of the energy coincides with the functional ap- 
proach in a special variational principle only for some equilibrium cases. 
The linear functional expressing the work of forces is related to the theory 
of distributions; a decomposition theorem associated with displacements 
(as test functions whose supports are C°° compact manifolds) uniquely de- 
termines a canonical zero order form (separated form) with respect both to 
the test functions and the transverse derivatives of contact test functions 9 . 
As Newton's principle is useless when we do not have any constitutive 
equation for the expression of forces, the linear functional method is useless 
when we do not have any constitutive assumption for the virtual work func- 
tional. The choice of the simple material theory associated with the Cauchy 
stress tensor corresponds with a constitutive assumption on its virtual work 
functional. It is important to notice that constitutive equations for the free 
energy \ an d constitutive assumption for the virtual work functional may 
be incompatible 10 : for any virtual displacement £ of an isothermal medium, 
the variation — 5\ must be equal to the virtual work of internal forces Srint- 
The equilibrium state is then obtained by the existence of a solution mini- 
mizing the free energy. 

The equation of motion of a continuous medium is deduced from the 
d'Alembert-Lagrange principle of virtual works which is an extension of 
the principle in mechanics of systems with a finite number of degrees of 
freedom: The motion is such that for any virtual displacement the virtual 
work of forces is equal to the virtual work of mass accelerations. 
Let us note: if the virtual work of forces is expressed in classical notations 
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in the form 

5t = J J J {f.C + tr[(-p l + 2/i VV).V<]}dw + J J T.C,ds (1) 

from the d'Alembert-Lagrange principle, we obtain not only the equations 
of balance momentum for a viscous fluid in the domain D but also the 
boundary conditions on the border S of D. We notice that expression (1) 
is not the Frechet derivative of any functional expression. 
If the free energy depends on the strain tensor F, then St must depend on 
VC and leads to the existence of the Cauchy stress tensor. If the free energy 
depends on the strain tensor F and on the overstrain tensor VF, then St 
must depend on V£ and V 2 £. 

Conjugated ( or transposed) mappings being denoted by asterisk, for any vec- 
tors a, b, we write a* b for their scalar product (the line vector is multiplied 
by the column vector) and ab* or a®b for their tensor product (the column 
vector is multiplied by the line vector). The product of a mapping A by a 
vector a is denoted by A a. Notation b* A means the covector c* defined 
by the rule c* = (A* b)*. The divergence of a linear transformation A is the 
covector divA such that, for any constant vector a, (div A) a = div (A a). 
We introduce a Galilean or fixed system of coordinates (x 1 , x 2 , x 3 ) which is 

also denoted by x as Euler or spatial variables. If / is a real function of x, 

df df df 

— — is the linear form associated with the gradient of / and — -r = (— ), ; 
ax ox 1 ax 

df 

consequently, ig-)* = g ra d/- The identity tensor is denoted by 1. 
Now, we present the method and its consequences in different cases of gradi- 
ent theory. As examples, we revisit the case of Laplace theory of capillarity 
and the case of van der Waals fluids. 



2. Virtual work of continuous medium 

The motion of a continuous medium is classically represented by a con- 
tinuous transformation ip of a three-dimensional space into the physical 
space. In order to describe the transformation analytically, the variables 
X =(X , X 2 , X 3 ) which single out individual particles correspond to mate- 
rial or Lagrange variables. Then, the transformation representing the mo- 
tion of a continuous medium is 

x = (f (X,t) or x i = ^(X 1 , X 2 , X 3 ,t), i = 1,2,3 

where t denotes the time. At t fixed the transformation possesses an inverse 
and continuous derivatives up to the second order except at singular sur- 
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faces, curves or points. Then, the diffcomorphism ip from the set Dq of the 
particles into the physical space D is an element of a functional space p of 
the positions of the continuous medium considered as a manifold with an 
infinite number of dimensions. 

To formulate the d'Alcmbcrt-Lagrange principle of virtual works, we in- 
troduce the notion of virtual displacements. This is obtained by letting 
the displacements arise from variations in the paths of the particles. Let a 
one-parameter family of varied paths or virtual motions denoted by {p>n\ 
and possessing continuous derivatives up to the second order and expressed 
analytically by the transformation 



with rj G O, where O is an open real set containing and such that 
$ (X,i; 0) = ip (X,i) or p = p (the real motion of the continuous medium 
is obtained when 77 = 0). The derivation with respect to 77 when 77 = is 
denoted by 5. Derivation 5 is named variation and the virtual displacement 
is the variation of the position of the medium 11 . The virtual displacement 
is a tangent vector to p in p (Sp G T v (p)). In the physical space, the 
virtual displacement Sp is determined by the variation of each particle: the 
virtual displacement of the particle x is such that £ = Sx when <5X = 0, 
Sn = 1 at T] = 0; we associate the field of tangent vectors to D 



where T X (D) is the tangent vector bundle to D at x. The concept of virtual 



Fig. 1. The boundary S of D is represented by a thick curve and its variation by a thin 
curve. Variation Sip of family {1^17} of varied paths belongs to T v (p), tangent space to 
(p) at ip. 



*(X,t;»?) 



X ED^C = iKx) = -T- |„= G T X (D) 
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work is purposed in the form: 

The virtual work is a linear functional value of the virtual displacement, 



where < . , . > denotes the inner product of 9 and Sip; then, 3 belongs to 
the cotangent space of p at <p (5 G T*{p)). 

In Relation (2), the medium in position p is submitted to the covector $> 
denoting all the stresses; in the case of motion, we must add the incrtial 
forces associated with the acceleration quantities to the volume forces. 
The d'Alembert-Lagrange principle of virtual works is expressed as: 



Theorem: // expression (2) is a distribution in a separated form, the 
d'Alembert-Lagrange principle yields the equations of motions and boundary 
conditions in the form 9 = 0. 

3. Some examples of linear functional of forces 

Among all possible choices of linear functional of virtual displacements, we 
classify the following ones: 

3.1. Model of zero gradient 
3.1.1. Model A.O 

The medium fills an open set D of the physical space and the linear func- 
tional is in the form 



where Fi (i = 1,2,3) denote the covariant components of the volume force 
F (including the inertial force terms) presented as a covector. The equation 
of the motion is 



St =< 3, 6(p > 



(2) 



For all virtual displacements, the virtual work is null. 

Consequently, representation (2) leads to: 

VS^e T V { P ), St = 




VxeD, F,=0^F = 



(3) 



3.1.2. Model B.O 



The medium fills a set D and the surface 5* is the boundary of D belonging 
to the medium; with the same notations as in section 3.1.1, the linear 
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functional is in the form 

Sr = J J J FiCdv + J f TiCda (4) 

Tj are the components of the surface forces (tension) T. From Eq. (4) , we 
obtain the equation of motion as in Eq. (3) and the boundary condition, 

VxeS, = T = 

3.2. Model of first gradient 

3.2.1. Model A.l 

With the previous notations, the linear functional is in the form 

*-///><*-<«)* 

where a\ = 1,2,3) are the components of the stress tensor a. Stokes 
formula gets back to the model B.O in the separated form 

where rij (j = 1, 2, 3) are the components of a covector which is the annu- 
lator of the vectors belonging to the tangent plane at the boundary S. It is 
not necessary to have a metric in the physical space; nevertheless, for the 
sake of simplicity it is convenient to use the Euclidian metric; the vector n 
of components n J (j = 1,2,3) represents the external normal to S relatively 
to D; the covector n* is associated with the components nj. We deduce 
the equation of motion 

VxeD, Fi + aij =0 <^ F + diver = (5) 

and the boundary condition 

VxeS, n 3 -aj = n* a = 

3.2.2. Model B.1/0: (Mixed model with first gradient in D and zero 
gradient on S) 

The linear functional is expressed in the form 

St = J J J (F i C i -oiC i j)dv + J JjiCds 
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Stokes formula yields the separated form 

Sr = J J J D { F *+ <o) Cdv + J J ft njaj) f ds (S. 0) 

and we deduce the equation of motion in the same form as Eq. ( 5) and the 
boundary condition 

VxeS, nj a\=Ti <S> n* a = T 

Model B.l / is the classical theory for elastic media and fluids in contin- 
uum mechanics. 

3.2.3. Model B.l 

The linear functional is expressed in the form 

5t = j j j D C - 4 C d ) dv + J J (Tt C + ll 0) ds (6) 

where the tensor 7 of components 7^ is a new term. The boundary of D is 
a surface S shared in a partition of TV parts S p of class C 2 , (p = 1,...,N) 
(Fig. 2). We denote by (i? m ) _1 the mean curvature of S; the edge T p of S p 
is the union of the limit edges T pq between surfaces S p and S q assumed to 
be of class C 2 and t is the tangent vector to T p oriented by n; n' is the unit 
external normal vector to T p in the tangent plane to S p : n'= t x n. Let us 
notice that: 

^0 = -7^C + ^ (7) 
where — -ff ( l ; consequently, from integration of the divergence of 




Fig. 2. The set D has a surface boundary S divided in several parts. The edge of S is 
denoted by T which is also divided in several parts with end points A m . 
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vector V on surfaces S p we obtain, 



Vids = -J J rij - Vjn l ^j ds + j n^di (8) 



We emphasize with the fact that Vj n l corresponds to the normal derivative 
dV 3 

to S„ denoted . An integration by parts of the term aj £*• in relation 

dn 

(6) and taking account of relations (7-8) implies 



St 



III F - v * + / I T - c ds+ IL L ^ d - + £ f Tr 



(S.I) 



with the following definitions 

Due to theorem 37 in 9 , the distribution (S .1) has a unique decomposition in 
displacements and transverse derivatives of displacements on the manifolds 
associated with D and its boundaries: expression (S.l ) is in a separated 
form. Consequently, the equation of motion is 

Vxefl, Ft = & F 1 = 

and the boundary conditions are 

V x G S, T} = 0, U = & T 1 = 0, L = 
Vx £ Tpq, Rpi+R q i = <=>■ R p + Rg = 

dC 

Term L is not reducible to a force: its virtual work Li—— is not the product 

dn 

of a force with the displacement the term L is an embedding action. 

3.3. Model of second gradient 
3.3.1. Mode^ A. 2 

The linear functional is in the form 



< ! J J D (F i C-oi(i j +Si k ( i tjk )dr 



Tensor S with Sj k = is an overstress tensor. An integration by parts 
of the last term brings back to the model B.l, 

5T = J J L ( Fl C (°* + ^ C ) dV + I L ^ C ' j dS 
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and the virtual work gets the separated form (S.l) with: 

Fl =Fi + cjL + S{ k jk volume force 

Tl = ~nj [a{ + S? fc " ^ (n k St) + ) surface force 

R P i = n'- n k Sj line force 

Li = rij n k Sf embedding action 

and consequently yields the same equation of motion and boundary condi- 
tions as in case B.l . 



3.3.2. Model B. 2 

The linear functional is in the form 

6T = J J L ^ + ^ fe ) dV+ I I ( TlC + ^ C ^ + U ^ k ) ds 

This functional yields two integrations successively on S p and on T pq with 
terms at the points A m . With obvious notations, for the same reasons as 
in section 3.2.3, the virtual work gets the separated form 

+ E / ( R p* ? + M Pi%) di + E <L is.2) 

p ^r p \ / m 

where C\. m (* = 1,2,3) are the components of C &t point A m . The calcula- 
tions are not expended. They introduce the curvature tensor on S p and the 
geodesic curvature of T pq _. Consequently, Fj , Tj, R pil (j) mi are associated 
with volume, surface, line and forces at points; Lj, Lf, M p i are embedding 
efforts of order 1 and 2 on S and of order 1 on the edge T. The equa- 
tion of motion and boundary conditions express that these seven tensorial 
quantities arc null on their domains of values D, S, T p and A m . 



4. Conclusion 

It is possible to extend the previous presentation by means of more complex 
medium with gradient of order n. The models introduce embedding effects 
of more important order on surfaces, edges and points. The (A. n) model 
refers to a (B.n-1) model: the fact that boundary surface S is (or is not) a 
material surface has now a physical meaning. Consequently, we can resume 
the previous presentation as follows: 
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a) The choice of a model corresponds to specify the part G of the algebraic 
dual T*(p) in which the efforts are considered: 3 e G C T*(p). 

b) In order to operate with the principle of virtual works and to obtain 
the mechanical equations in the form 9 = 0, it is no matter that the part 
G of the dual is separating (V3 6 G, < 5, 5<p > = =>■ 5(p = 0), but it is 
important the part G is separated (Sg G,M 5<p G T v (p), < 3, Sip > = => 
3 = 0). 

c) The functionals A.l , B.l / 0, A. 2, B.2 are not separated: if 3 consists in 
the data of the fields F, a, T, it is not possible to conclude that the fields 
are zero. 

d) Functionals in A.0, B.0, S.l, S.2 ... are separated: if the fields S 1 , T 1 , 
R. , L, . . . are continuous then, by using the fundamental lemma of variation 
calculus, their values must be equal to zero. They are the only functionals 
we must know for using the principle of virtual works; it is exactly as for a 
solid: the torque of forces is only known in the equations of motion. 

e) When the fields are not continuous on surfaces or curves, we have to 
consider a model of greater order in gradients and to introduce integrals on 
inner boundaries of the medium. 

For conservative medium, the first gradient theory corresponds to the 
compressible case. The theory of fluid, elastic, viscous and plastic media 
refers to the model (S.0). The Laplace theory of capillarity in fluids refers 
to the model (S.l). To take into account superficial effects acting between 
solids and fluids, we use the model of fluids endowed with capillarity (S.2); 
the theory interprets the capillarity in a continuous way and contains the 
Laplace theory of capillarity; for solids, the model corresponds to "elastic 
materials with couple stresses" indicated by Toupin in 12 . 

5. Example 1: The Laplace theory of capillarity 

Liquid- vapor and two-phase interfaces are represented by a material surface 
endowed with an energy relating to Laplace surface tension. The interface 
appears as a surface separating two media with its own characteristic be- 
havior and energy properties 13 (when working far from critical conditions, 
the capillary layer has a thickness equivalent to a few molecular beams 14 ). 
The Laplace theory of capillarity refers to the model B.l in the form (S.l) 
as following: for a compressible fluid with a capillary effect on the wall 
boundaries, the free energy is in the form 
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where a(p) is the fluid specific energy, p is the matter density and coeffi- 
cients a p are the surface tensions of each surface S p . Surface integrations 
are associated to the space metric; the virtual work of internal forces is 

ST int = J J J D -P,iCdv+J2j j s n t C>**+J2j T a v<Cdt 

where p = p 2 a'{p) is the fluid pressure. The external force (including inertial 
forces) is the body force pi defined in D, the surface force is T defined on S 
and the line force is R defined on F D'Alembert-Lagrange principle yields 
the equation of motion and boundary conditions: 

Vx e D, -p ti + pfi = -gradp + pf = 0, 

ft ■ n 

VxeS, p n, + T, - a p — - = <^ pn + T - a p — = 

VxeTpq, OpTi'pi + a g n' gi + Ri = ^> a p n' p + a q n' q + R = 
Boundary conditions are Laplace equation and Young-Dupre condition. 

6. Example 2: Fluids endowed with internal capillarity 

For interfacial layers, kinetic theory of gas leads to laws of state associated 
with non-convex internal energies 15,16 . This approach dates back to van 
der Waals 17 , Korteweg 18 , corresponds to the Landau-Ginzburg theory 19 
and presents two disadvantages. First, between phases, the pressure may 
become negative; simple physical experiments can be used to cause traction 
that leads to these negative pressure values 20 . Second, in the field between 
bulks, internal energy cannot be represented by a convex surface associated 
with density and entropy; this fact seems to contradict the existence of 
equilibrium states; it is possible to eliminate this disadvantage by writing 
in an anisotropic form the stress tensor of the capillary layer which allows 
to study interfaces of non- molecular size near a critical point. 
One of the problems that complicates this study of phase transformation 
dynamics is the apparent contradiction between Korteweg classical stress 
theory and the Clausius-Duhem inequality 21 . Proposal made by Eglit 22 , 
Dunn and Serrin 23 , Casal and Gouin 24 and others rectifies this anomaly 
for liquid- vapor interfaces. The simplest model in continuum mechanics 
considers a free energy as the sum of two terms: a first one correspond- 
ing to a medium with a uniform composition equal to the local one and a 
second one associated with the non-uniformity of the fluid 15,17 . The sec- 
ond term is approximated by a gradient expansion, typically truncated to 
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the second order. The model is simpler than models associated with the 
renormalization-group theory 25 but has the advantage of easily extending 
well-known results for equilibrium cases to the dynamics of interfaces 26-28 . 
We consider a fluid D in contact with a wall 5. Physical experiments prove 
that the fluid is nonhomogeneous in the neighborhood of S 16 . The internal 
energy e is also a function of the entropy. In the case of isothermal motions, 
the internal energy is replaced by the free energy. In the mechanical case, 
the entropy and the temperature are not concerned by the virtual displace- 
ments of the medium. Consequently, for isentropic or isothermal motions, 
e = f(p,0) where = (gradp) 2 . The fluid is submitted to external forces 
represented by a potential O as a function of Eulerian variables x. To obtain 
boundary conditions it is necessary to know the wall effect. An explicit form 
for the energy of interaction between surfaces and liquids is proposed in 29 . 
We denote by B the surface density of energy at the wall. The total energy 
E of the fluid is the sum of three potential energies: Ef (bulk energy), E p 
(external energy) and E$ (surface energy). 

E f = J J J D Pe(p,0) dv, E p = J J J P^) dv, E s = J J^B ds 
We have the results (see Appendix): SEf = J J J (— div a) C dv 

+ + G£ + grad; »- 4 + nV ) c ) ds - 1 A c M 

with cr = -Pl-Cgradp®gradp= -Pl-C (|^)* where C = 2pe' n , 

ax ax p 

P = p 2 e' p — p div(C gradp), e' p (or s'g) denoting the partial derivative of e 

dp dp dp 

with respect to p (or 0), ( n = n* A = Cp — where — = — — n and 

an an ax 

grad tff denotes the tangential part of the gradient relatively to S. 



5E p = J J J p dv = J J J p (grad* fi) C dv; and, 

SE S = J J s { SB - (jf n * + S^d* tg B^j C^ds + J^B n'* C dJt 

The density in the fluid has a limit value p s at the wall S and B is assumed 
to be a function of p s only 29 . Then, SB = B'(p s )5p s = —p s B'(p s )dWC, 
where div£ is computed on S 11 . Let us denote G = —p s B'(p s ); Appendix 
yields 



J J 5Bds = j j GdWCds 
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, , G t-(S n " +g " d « G )4 dS + /r Gn '"^ f 

d(n _ f2H_ 
is I dn \R m / j Jr 

with H = B(p 3 )- p s B'(p s ). Then, 

= (pgrad*ft-div<7)C<fo - J (A- H)n'* Cd£ 

+ IL {G -" A) t + (^4^ »* + ^ - h > + -'") < «» <9) 

At equilibrium, <5r = — <5_E = 0. The fundamental lemma of variation cal- 
culus associated with separated form (9) corresponding to (S.2), yields: 

6.1. Equation of equilibrium 

From any arbitrary variation xeJ)-t((x) such that ( = on 5, we get 
J J J (P grad*Q - diver) £ ds = 0. Then, 

pgrad*f2- diver = (10) 

This equation is written in the classical form of equation of equilibrium 24 . 
It is not the same for the boundary conditions. 

6.2. Boundary conditions 
6.2.1. Case of a rigid wall 

We consider a rigid wall; on S, the virtual displacements satisfy the condi- 
tion n* C = 0. Then, at the rigid wall, V x e S -> £(x) such that n* C = 0, 



(G - ^4)-^7 + ( -5-5 ^ n* + grad t ;(A -ff) + nV]Crfs = 



s 

Due to fj = cr*, we deduce the boundary conditions (11-12) 

VxeS, G-A = 0, (11) 
and there exists a Lagrange multiplier x e 5 — > A(x) e R such that, 
2(A - H) 

The edge T of S belongs to the solid wall and consequently on T, £ = ryt: 
the integral on T is null and does not yield any additive condition. 



VxeS 1 , -i-g ^ n + grad tg (A -i?) + <7n = An (12) 
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6.2.2. Case of an elastic wall 

The equilibrium equation (10) is unchanged. On S, the condition (11) is 
also unchanged. The only different condition comes from the fact that we 
do not have anymore the slipping condition for the virtual displacement 
on S, (n* £ = 0). Due to the possible deformation of the wall, the virtual 

work of stresses on S is SE e = J J k*£ ds + J H*£d£ where k = a e n 

is the stress (loading) vector associated with stress tensor o~ e of the elastic 
wall and R is the line force due to the elasticity of the line. Relation (12) 
is replaced by 



VxeS, 2 ^— ^ '- n + grad 4g (A -H)+an + K = Q 



(A_H) 

We obtain an additive condition on T in the form (H — A) n' + R = and 
due to condition (11), 



VxeT, Bn'+R=0 (13) 
(If T is the union of edges T p , B n' is replaced by ^ B p n' p ). 



6.3. Analysis of the boundary conditions 

Eq. (11) yields C ^ + B'{p s ) = 0; the definition of a implies: 

a n = -Pn — C ^ grad p 
dn 

Due to the fact that the tangential part of Eq. (12) is always verified, 
the only condition comes from Eq. (11); Eq. (12) yields the value of the 
Lagrange multiplier A and Eq. (13) the value of R. For an clastic (non- 
rigid) wall we obtain, 

K tg = and n n =P+^--B'(p s )^ (14) 

Km dn 



where n tg and K )i RFC the tangential and the normal components of k. 
Taking into account of Eq. (14) wc obtain the stress values at the non rigid 
elastic wall. The surface energy is 29 : B(p s ) — —71 p s + — p s where 71 and 
72 are two positive constants and the fluid density condition at the wall is 



n d P 
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If we denote by p B = 71 / 72 the bifurcation fluid density at the wall, due to 
the fact C is positive constant 16 , we obtain: if p s < p B , (or p s > p B ), 

is positive ( or negative) and we have a lack ( or excess) of fluid density at 
the wall. Such media allow to study fluid interfaces and interfacial layers 
between fluids and solids and lead to numerical and asymptotic methods 30 . 
The extension to the dynamic case is straightforward: Eq. (10) yields 

pT* — diva + p grad*f2 = 

Vector r is the acceleration; boundary conditions (11-14) are unchanged. 
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Appendix 

Let S be a surface in the 3-dimensional space and n its external normal 
extended locally in the vicinity of S by the expression n(x) = grad g?(x), 
where d is the distance of a point x to S; for any vector field w, we obtain 31 : 

/nt ,. c* n dw 

rotln x w) = n div w — w div n + — — w — — n 

ax ax 

du 2 

From n*— — = and divn = wc deduce on S, 

ax R m 

2 5w 
n*rot(n x w) = div w H n*w n* — — n (15) 

Rm OX 



We deduce: for any scalar field A and w = A C 

dpa 2A 
dn R 



AdivC = A^-—Cn- (grad t * s A) < + n*rot (An x C) (16) 



Let us calculate 5Ef, D is a material volume, then 5Ef = J J J p Se dv 

>4''-^'£(£)'- c (£-si)(iy 

= div(C gradp 5p) — div(C gradp) Sp — tr yCgr&dp grad*p ^ j 
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Due to Sp = — pdivC (see ) 



SE f = J J J + div(C grad p grad*p)^ Cdv 

III ^ V ^ /9 ^ ra< ^ /9 div£ + Cgradp grad* p C + PC) dv 
= j J J -(diver) Cdv + J J (-AdivC + nVC)ds 

From Eq. (16), we deduce immediatly: SEf = J J J (— divrj)£ dv + 

~ A ^ + {it a ' + A+ nV ) <} ds ~ I An " <* 



Let us calculate SEg; due to E$ = J J B dct (n, rfix, <i 2 x) where 
(iix and g? 2 x are two coordinate lines of S, we get: 

E s = I I B det F det (-F _1 n, rfiX, ri 2 X) 

J JSa 

where Sq is the image of S in a reference space with Lagrangian coordinates 

(9x f c^x* 1 

X and F is the deformation gradient tensor 7^- of components < r • 

Then, 5E S = [ [ SB dct F dct (F _1 n, ciiX, d 2 X) 

J JSa 

+ [ [ BS (det F det(F- 1 n,d 1 X,d 2 X)) . 

J JSa 

with / / BS (det F det (F- 1 n,d 1 X,d 2 X)) = 

J JSn 

(dn 
— C, dix, d 2 x 

-Bdct ^n,dix,rf 2 xj = J j (dxv(BQ - (grad*B)C-Bn*^nJ ds 

Relation (15) yields: div(BC) + -^- n *C - n *^T ln = n * rot ( B n x C), 
/ / BS (dct Fdct(F- 1 n,d 1 X,d 2 X)) = 
— n* + grad*5 (nn* -l))£ds+ / n* rot (J3nx() ds 



So 



R 



where grad*B (nn* — 1) belongs to the cotangent plane to S; we obtain 



2B 
FL 



SE S = j j (SB- [ ^- n*+ grad* ff B) C ) ds + \Bn'\d£. 
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